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The real-time dynamics of a classical spin in an external magnetic field and locally exchange 
coupled to an extended one-dimensional system of non-interacting conduction electrons is studied 
numerically. Retardation effects in the coupled electron-spin dynamics are shown to be the source 
for the relaxation of the spin in the magnetic field. Total energy and spin is conserved in the 
non-adiabatic process. Approaching the new local ground state is therefore accompanied by the 
emission of dispersive wave packets of excitations carrying energy and spin and propagating through 
the lattice with Fermi velocity. While the spin dynamics in the regime of strong exchange coupling 
J is rather complex and governed by an emergent new time scale, the motion of the spin for 
weak J is regular and qualitatively well described by the Landau-Lifschitz-Gilbert (LLG) equation. 
Quantitatively, however, the full quantum-classical hybrid dynamics differs from the LLG approach. 
This is understood as a breakdown of weak-coupling perturbation theory in J in the course of time. 
Furthermore, it is shown that the concept of the Gilbert damping parameter is ill-defined for the 
case of a one-dimensional system. 

PACS numbers: 75.78.-n, 75.78.Jp, 75.60.Jk, 75.10.Hk, 75.10.Lp 


I. INTRODUCTION 

The Landau-Lifshitz-Gilbert (LLG) equatiorP^ has 
originally been considered to describe the dynamics of 
the magnetization of a macroscopic sample. Nowadays it 
is frequently used to simulate the dynamics of many mag¬ 
netic units coupled by exchange or magnetostatic interac¬ 
tions, i.e., in numerical micromagnetics.l^The same LLG 
equation can be used on an atomistic level as well.^SHU Por 
a suitable choice of units and for several spins Sm{t) at 
lattice sites m, it has the following structure: 

= S^t) XB + J2 X S„it) 

n 

+ ( 1 ) 

n 

It consists of precession terms coupling the spin at site 
m to an external magnetic field B and, via exchange 
couplings Jmm to the spins at sites n. Those pre¬ 
cession terms typically have a clear atomistic origin, 
such as the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
interactiorPsHUl which is mediated by the magnetic po¬ 
larization of conduction electrons. The non-local RKKY 
couplings Jmn = J'^Xmn are given in terms of the ele¬ 
ments Xmn of the static conduct ion-electron spin suscep¬ 
tibility and the local exchange J between the spins and 
the local magnetic moments of the conduction electrons. 
Other possibilities comprise direct (Heisenberg) exchange 
interactions, intra-atomic (Hund’s) couplings as well as 
the spin-orbit and other anisotropic interactions. The re¬ 
laxation term, on the other hand, is often assumed as lo¬ 
cal, amn = ^mnOi, and represented by purely phenomeno¬ 
logical Gilbert damping constant a only. It describes the 
angular-momentum transfer between the spins and a usu¬ 
ally unspecified heat bath. 


On the atomistic level, the Gilbert damping must be 
seen as originating from microscopic couplings of the 
spins to the conduction-electron system (as well as to 
lattice degrees of freedom which, however, will not be 
considered here). There are numerous studies where the 
damping constant, or tensor, a has been computed nu¬ 
merically from a more fundamental model including elec¬ 
tron degrees of freedom explicitljff^®^ or even from first 
principles.liSHSI] All these studies rely on two, partially 
related, assumptions: (i) The spin-electron coupling J 
is weak and can be treated perturbatively to lowest or¬ 
der, i.e., the Kubo formula or linear-response theory is 
employed, (ii) The classical spin dynamics is slow as 
compared to the electron dynamics. These assumptions 
appear as well justified but they are also necessary to 
achieve a simple effective spin-only theory by eliminat¬ 
ing the fast electron degrees of freedom. 

The purpose of the present paper is to explore the 
physics beyond the two assumptions (i) and (ii). Us¬ 
ing a computationally efficient formulation in terms of 
the electronic one-particle reduced density matrix, we 
have set up a scheme by which the dynamics of classi¬ 
cal spins coupled to a system of conduction electrons can 
be treated numerically exactly. The theory applies to ar¬ 
bitrary coupling strengths and does not assume a separa¬ 
tion of electron and spin time scales. Our approach is a 
quantum-classical hybrid theorj^^ which may be charac¬ 
terized as Ehrenfest dynamics, similar to exact numerical 
treatments of the dynamics of nuclei, treated as classical 
objects, coupled to a quantum system of electrons (see, 
e.g.. Ref. [23l for an overview). Some other instructive ex¬ 
amples of quantu m-cla ssical hybrid dynamics have been 
discussed recently.l^H^ 

The obvious numerical advantage of an effective spin- 
only theory, as given by LLG equations of the form Q, 
is that in solving the equations of motion there is only 
the time scale of the spins that must be taken care of. As 






2 


compared to our hybrid theory, much larger time steps 
and much longer propagation times ca n be achieved. Op¬ 
posed to ab-initio approache d^^ l ^^ * ^^ we therefore con¬ 
sider a simple one-dimensional non-interacting tight- 
binding model for the conduction-electron degrees of free¬ 
dom, i.e., electrons are hopping between the nearest- 
neighboring sites of a lattice. Within this model ap¬ 
proach, systems consisting of about 1000 sites can be 
treated easily, and we can access sufficiently long time 
scales to study the spin relaxation. An equilibrium state 
with a half-filled conduction band is assumed as the ini¬ 
tial state. The subsequent dynamics is initiated by a 
sudden switch of a magnetic field coupled to the classical 
spin. The present study is performed for a single spin, 
i.e., we consider a classical-spin Kondo-impurity model 
with antiferromagnetic local exchange coupling J, while 
the theory itself is general and can be applied to more 
than a single or even to a large number of spins as well. 

As compar ed t o the conventional (quantum-spin) 
Kondo model,^23111 model considered here does not 
account for the Kondo effect and therefore applies to sit¬ 
uations where this is absent or less important, such as 
for systems with large spin quantum numbers S, strongly 
anisotropic systems or, as considered here, systems in a 
strong magnetic field. To estimate the quality of the 
classical-spin approximation a priori is difficult Por 
one-dimensional systems, however, a quantitative study 
is possible by comparing with full quantum calculations 
and will be discussed elsewhere. 

There are different questions to be addressed: For 
dimensional reasons, one should expect that linear- 
response theory, even for weak J, must break down at 
long times. It will therefore be interesting to compare 
the exact spin dynamics with the predictions of the LLG 
equation for different J. Furthermore, the spin dynamics 
in the long-time limit can be expected to be sensitively 
dependent on the low-energy electronic structure. We 
will show that this has important consequences for the 
computation of the damping constant a and that a is 
even ill-defined in some cases. An advantage of a full 
theory of spin and electron dynamics is that a precise 
microscopic picture of the electron dynamics is available 
and can be used to discuss the precession and relaxation 
dynamics of the spin from another, namely from the elec¬ 
tronic perspective. This information is in principle exper¬ 
imentally accessible to spin-resolved scanning-tunnelling 
microscope techniqued^SHsej important for an at om- 
istic understanding of nano-spintronics devicesWe 
are particularly interested in the physics of the system 
in the strong-J regime or for a strong field B where the 
time scales of the spin and the electron dynamics become 
comparable. This has not yet been explored but could 
become relevant to understand real-time dynamics in re¬ 
alizations of strong-J Kondo-lattice models by means of 
ultraco ld ferm ionic Yb quantum gases trapped in optical 
lattices.l^^*^ 

The paper is organized as follows: We first introduce 
the model and the equations of motion for the exact 


quantum-classical hybrid dynamics in Sec. [ll|an d discuss 
some computational details in Sec. |III| Sec. |IV| provides 
a comprehensive discussion of the relaxation of the clas¬ 
sical spin after a sudden switch of a magnetic field. The 
reversal time as a function of the interaction and the field 
strength is analyzed in detail. We then set the focus on 
the conduction-electron system which induces the relax¬ 
ation of the classical spin by dissipation of energy. In Sec. 
[Vl the linear-response approach to integrate out the elec¬ 
tron degrees of freedom is carefully examined, including a 
discussion of the additional approximations that are nec¬ 
essary to re-derive the LLG equation and the damping 
term in particular. Sec. VI summarizes the results and 
the main conclusions. 


II. MODEL AND THEORY 

We consider a classical spin S with [S'! = 1/2, which is 
coupled via a local exchange interaction of strength J to 
the local quantum spin at the site iq of a system of N 
itinerant and non-interacting conduction electrons. The 
conduction electrons hop with amplitude —T (T > 0) 
between non-degenerate orbitals on nearest-neighboring 
sites of a D-dimensional lattice, see Fig. L is the 
number of lattice sites, and n = N/L is the average 
conduction-electron density. 

The dynamics of this quantum-classical hybrid 
systeirP^ is determined by the Hamiltonian 

H = -TY, + Js.,S - BS . (2) 

Here, Cia annihilates an electron at site i = 1, ...,L with 
spin projection a =t,4., and Si = \ is 

the local conduction-electron spin at z, where cr denotes 
the vector of Pauli matrices. The sum runs over the 
different ordered pairs {ij) of nearest neighbors. B is 
an external magnetic field which couples to the classical 
spin. 

To be definite, an antiferromagnetic exchange coupling 
J > 0 is assumed. If S was a quantum spin with 
S = 1/2, E g. would represent the single-impurity 
Kondo model.123^ However, in the case of a classical spin 
considered here, there is no Kondo effect. The semiclas- 
sical single-impurity Kondo model thus applies to sys¬ 
tems where a local spin is coupled to electronic degrees 
of freedom but where the Kondo effect absent or sup¬ 
pressed. This comprises the case of large spin quantum 
numbers S, or the case of temperatures well above the 
Kondo scale, or systems with a ferromagnetic Kondo cou¬ 
pling J < 0 where, for a classical spin, we expect a qual¬ 
itatively similar dynamics as for J > 0. 

We assume that initially, at time t = 0, the clas¬ 
sical spin S(t = 0) has a certain direction and that 
the conduction-electron system is in the corresponding 
ground state, i.e., the conduction electrons occupy the 
lowest N one-particle eigenstates of the non-interacting 
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Hamiltonian Eq. ([^ for the given S' = = 0) up to 

the chemical potential A non-trivial time evolution is 
initiated if the initial direction of the classical spin and 
the direction of the field B are non-collinear. 

To determine the real-time dynamics of the electronic 
subsystem, it is convenient to introduce the reduced one- 
particle density matrix. Its elements are defined as ex¬ 
pectation values, 

Pii\aa '= {^i'a'^ ’ (^) 

in the system’s state at time t. At t = 0 we have p(0) = 
0(/r — T{0)). The elements of p(0) are given by 

PiiT,i'cr' (0) = y ] Uia-^k^ip ~ ^k)U^k,i'a' ’ ('^) 

k 

where 0 is the step function and where U is the uni¬ 
tary matrix diagonalizing the hopping matrix T(0), i.e., 
U^T{0)U = e with the diagonal matrix e given by the 
eigenvalues of T(0). The hopping matrix at time t is 
can be read off from Eq. (H- It comprises the physical 
hopping and the contribution resulting from the coupling 
term. Its elements are given by 

T era' • (b) 

Here = 1 if i, z' are nearest neighbors and zero else. 

There is a closed system of equations of motion for the 
classical spin vector S {t) and for the one-particle density 
matrix p(t). The time evolution of the classical spin is 
determined via {d/dt)S{t) = {S, } by the classical 

Hamilton function = (H). This equation of mo¬ 

tion is the only known way to consistently describe the 
dynamics of quantum-classical hybrids (see Refs. l22l41l42l 
and references therein for a general discussion). The 
Poisson bracket between arbitra ry fun ctions A and B of 
the spin components is given by,S2Hll 

r A 7-jT ry //^\ 

> (6) 

where the sums run over x, y, z and where Eap-y is the 
fully antisymmetric e-tensor. With this we find 

j^Sit) = J{s,„)txSit)-BxSit). ( 7 ) 

This is the Landau-Lifschitz equation where the expec¬ 
tation value of the conduction-electron spin at ig is given 

by 


{^io)t — 2 'y / PiorT,io<7’ jt) O'o-'g- , 
(Tcr' 


( 8 ) 


and where J{sig)t acts as an effective time-dependent 
internal field in addition to the external field B. 


B 



io 


FIG. 1: (Color online) Classical spin S{t) coupled via an an¬ 
tiferromagnetic local exchange interaction of strength J to a 
system of conduction electrons hopping with nearest-neighbor 
hopping amplitude T over the sites of a one-dimensional lat¬ 
tice with open boundaries. The spin couples to the central 
site io of the system and is subjected to a local magnetic field 
of strength B. 


The equation of motion for {si)t reads as 

n.n. ^ 

+ (9) 

j na' 

where the sum runs over the nearest neighbors of i. The 
second term on the right-hand side describes the coupling 
of the local conduct ion-electron spin to its environment 
and the dissipation of spin and energy into the bulk of 
the system (see below). Apparently, the system of equa¬ 
tions of motion can only be closed by considering the 
complete one-particle density matrix Eq. (§. It obeys a 
von Neumann equation of motion, 

z|p(t) = [T(t),p(t)] (10) 

as is easily derived, e.g., from the Heisenberg equation of 
motion for the annihilators and creators. 

As is obvious from the equations of motion, the real¬ 
time dynamics of the quantum-classical Kondo-impurity 
model on a lattice with a finite but large number of sites L 
can be treated numerically exactly (see also below). Nev¬ 
ertheless, the model comprises highly non-trivial physics 
as the electron dynamics becomes effectively correlated 
due to the interaction with the classical spin. In addi¬ 
tion, the effective electron-electron interaction mediated 
by the classical spin is retarded: electrons scattered from 
the spin at time t will experience the effects of the spin 
torque exerted by electrons that have been scattered from 
the spin at earlier times t' <t. 


III. COMPUTATIONAL DETAILS 


Eqs. (§, 0, 0 and ( [l0| ) represent a coupled non¬ 
linear system of first-order ordinary differential equations 
which can be solved numerically. By blocking up the 
von Neumann equation (10), the differential equations 
are written in a standard form y = f{y{t),t)y where y{t) 
is a high-dimensional vector, such that an explicit Runge- 
Kutta method can be applied. A high-order propagation 
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technique is usecP^ which provides the numerically exact 
solution up to 6-th order in the time step At. For a typ¬ 
ical system consisting of about L = 10^ sites this implies 
that ~ 10® coupled equations are solved. 

We consider a one-dimensional system with open 
boundaries consisting of L = 1001 sites and a local per¬ 
turbation at the central site *o of the system, see Fig. 
For a half-filled tight-binding conduction band the Fermi 
velocity up = 2T roughly determines the maxi mum s peed 
of the excitations and defines a “light cone” .ESMI This 
means that finite-size effects due to scattering at the sys¬ 
tem boundaries become relevant after a propagation time 
imax 500 (in units of 1/T). A time step At = 0.1 is 
usually sufficient for reliable numerical results up to tmaxj 
i.e., about 5000 time steps are performed. The compu¬ 
tational cost is moderate, and calculations can be per¬ 
formed in a few hours on a standard desktop computer. 

Assuming, for example, that B = (0,0, B), the Hamil¬ 
tonian is invariant under rotations around the z axis. It 
is then easily verified that not only the length of the spin 
[S'! = 1/2 is conserved but also the total number of con¬ 
duction electrons, 

Ntot = ) ( 11 ) 

i(7 

the z-component of the total spin, 

'S'tot.z = Sz + , (12) 

i 

as well as the total energy, 

Etot = {H) = tr (p(t)r(t)) - BSit) . (13) 

Despite the fact that the model does not include a di¬ 
rect (e.g.. Coulomb) interaction among the conduction 
electrons, the average occupation numbers of the basis of 
one-particle states in which the hopping matrix T{t) is 
diagonal at time t are not conserved. This is due to the 
effective retarded interaction mediated by the classical 
spin. Hence, the system is not integrable, unlike a free 
fermion gas. The conservation of the above-mentioned 
global observables serves as a sensitive check for the ac¬ 
curacy of the numerical procedure. 


IV. COUPLED SPIN AND ELECTRON 
DYNAMICS 

A. Spin relaxation 

Fig-i shows the real-time dynamics of the classical 
spin for J = 1. Energy and time units are fixed by the 
nearest-neighbor hopping T = 1 throughout the paper. 
Initially, for t = 0, the spin is oriented (almost) antipar¬ 
allel to the external local field B — (0,0, B) with B = 1, 
i.e., initially 5'x(0) = ^sind, Sy{0) = 0, 5'2(0) = —^cosi? 



FIG. 2: (Color online) Real-time dynamics of the classical 
spin. Upper panel: Bloch sphere representation. Lower panel: 
X and z component of S{t) (|S’| = 1/2). Calculations for 
exchange coupling J = 1 and field strength B = 1 and for 
a system of L = 1001 sites. (L — 1001 is kept fixed for the 
rest of the paper). Energy and time units are fixed by the 
nearest-neighbor hopping T = 1. 


where a non-zero but small polar angle = 7r/50 is nec¬ 
essary to slightly break the symmetry of the initial state 
and to start the dynamics. 

For the same setup, the Landau-Lifschitz-Gilbert equa¬ 
tion would essentially predict two effects: first, a preces¬ 
sion of the classical spin around the field direction with 
Larmor frequency ujl = B^, and second, a relaxation 
of the spin to the equilibrium state with S parallel to 
B = Bez for t —>■ oo. Both effects are also found in 
the full dynamics of the quantum-classical hybrid model. 
The frequency of the oscillation of Sx{t) that is seen in 
Fig.@ is wl , and Sz is reversed after a few hundred time 
units. The precessional motion is easily explained by the 
torque on the spin exerted by the field according to Eq. 
0. The explanation of the damping effect is more in¬ 
volved: 

Even for the high field strength considered here, the 
spin dynamics is slow as compared to the characteristic 
electronic time scale such that it could be reasonable to 
assume the electronic system being in its instantaneous 
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exchange coupling J 


FIG. 3: (Color online) Time dependence of the angle 'y{t) 
enclosed by S{t) and {sio)t for B = 0.1 and different J as 
indicated. 


FIG. 4: (Color online) Time for a spin reversal Sz = | cos( 7 r — 
1 ?) = — I cos(’d) —>■ Sz = I cos(i?) for = i9i = tt/SO (reversal 
time Ti) and for i9 = 1^2 = 7r/25 (reversal time T 2 ). Calcula¬ 
tions as a function of J for fixed B — 0.1. 


ground state at any instant of time and corresponding 
to the configuration of the classical spin. This, however, 
would imply that the expectation value {sig)t of the local 
conduction-electron moment at zq is always strictly par¬ 
allel to S{t) and, hence, there would not be any torque 
mediated by the exchange coupling J on S{t). 

In fact, the direction of {si^)t is always somewhat be¬ 
hind the “adiabatic direction”, i.e., behind —S{t): This 
is shown in Fig. |^for a field of strength B = 0.1 where 
the spin dynamics is by a factor 10 slower, compared to 
Fig.H and for different stronger exchange couplings J. 
Even in this case the process is by no means adiabatic, 
and the angle ^{t) between S{t) and {sig)t is close to but 
clearly smaller than 7 = tt at any instant of time. This 
non-adiabaticity results from the fact that the motion 
of the classical spin affects the conduction electrons in 
a retarded way, i.e., it takes a finite time until the local 
conduction-electron spin at ig reacts to the motion 

of the classical spin. 

This retardation effect results in a torque J{sig)t x 
S{t) 0 exerted on the classical S (t) in the -\-z direction 
which adds to the torque due to B and which drives the 
spin to its new equilibrium direction. Hence, retardation 
is the physical origin of the Gilbert spin damping. 

With increasing time, the deviation of 'y{t) from the 
instantaneous equilibrium value 7 = tt increases in mag¬ 
nitude, i.e., the direction of {sig)t is more and more be¬ 
hind the adiabatic direction, and the torque increases. Its 
magnitude is at a maximum at the same time when the 
oscillating x,y components of S{t) are at a maximum 
(see Fig. [^. The z-component of the torque does not 
vanish before the spin has reached its new equilibrium 
position S{t) (x e^. 

Fig.i shows results for different J. Generally, non- 
adiabatic effects show up if the typical time scale of the 
dynamics is faster than the relaxation time, i.e., the time 
necessary to transport the excitation away from the lo¬ 
cation ig where it is created initially. Roughly this time 
scale is set by the inverse hopping 1/T. One therefore 
expects that, for fixed T, a stronger J implies a stronger 


retardation of the conduction-electron dynamics. The re¬ 
sults for different J shown in Fig. in fact show that the 
maximum deviation of ■j(t) from the adiabatic direction 
7 = TT increases with increasing J (for very strong J the 
dynamics becomes much more complicated, see below). 
This results in a stronger torque on S{t) in z direction 
and thus in a stronger damping. The picture is also qual¬ 
itatively consistent with the LEG equation as the Gilbert 
damping constant a increases with J. 

The spin (almost) reverses its direction after a finite 
reversal time r which is shown in Fig. as a function 
of J. Calculations have been performed for an initial 
direction of the classical spin with 5'z(0) = —( 1 / 2 ) cos d 
with two different polar angles di_ 2 - If is sufficiently 
small, the results for different z? are expected to differ by 
a J and B independent constant factor only. As Fig. 
demonstrates, the ratio ti/t 2 is in fact nearly constant. 
For weak J and up to coupling strengths of about J < 30, 
we find that the reversal time decreases with increasing 
J. With increasing J, the retardation effect increases, as 
discussed above, and the stronger damping results in a 
shorter reversal time. 



FIG. 5: (Color online) Reversal time ti as function of B for 
fixed J — 2 and J = 4 as indicated. 
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The prediction of the LLG equation for the reversal 
time of a single spin r can be derived analyticallj^^ and 
is given by 


1 + aM 1/2-5,(0) 

a B 1/2+ 5,(0) 


(14) 


However, down to the smallest J for which r can be cal¬ 
culated reliably, our results for the full spin dynamics do 
not scale as t oc 1/J^ as one would expect for w eak J 
assuming that a oc. (see discussion in Sec. VB). 

The B dependence of the reversal time is shown in Fig. 

With increasing field strength, the classical spin S{t) 
precesses with a higher Larmor frequency ujl ~ B around 
the z axis. Hence, non-adiabatic effects increase. The 
stronger the field, the more delayed is the precessional 
motion of the local conduction-electrons spin {sig)t- This 
results in a stronger torque in +z direction exerted on 
the classical spin. Therefore, the relaxation is faster and 
the reversal time r smaller. For weak and intermediate 
field strengths, r is roughly proportional to 1/B. This is 
consistent with the prediction of the LLG equation, see 


Eq. (14). 


In the limit of very strong fields one would expect an 
increase of the reversal time with increasing B since the 
field term will eventually dominate the dynamics, i.e., 
only the precessional motion survives which implies a di¬ 
verging reversal time. In fact, for a field strength exceed¬ 
ing a critical strength B^,, which depends on J, there is 
no full relaxation any longer, and r = oo. This strong- 
B regime cannot be captured by the LLG equation and 
deserves further studies. 

The strong-J regime is interesting as well. For cou¬ 
pling strengths exceeding J « 30 the reversal time in¬ 
creases with J (see Fig. |^. Eventually, the reversal time 
must even diverge. This is obvious as the dynamics is 
described by a simple two-spin model in the limit J = oo 
which cannot show spin relaxation. The corresponding 
equations of motion are obtained by from Eqs. ([^ and 
([^ by setting t = 0 and B = 0: 


-Sit) = X S{t) , 

^(s*o)t = JSit) X 


(15) 


Note that we have |(sio)t| = 1/2 for J —>• oo. The equa¬ 
tions are easily solved by exploiting the conservation of 
the total spin Stot = Sl{t) (s^o)*. Both spins process 
with constant frequency luq = JStot around Stot- Their 
components parallel to 5tot are equal, and their compo¬ 
nents perpendicular to 5tot are anti-parallel and of equal 
length. 

However, the two-spin dynamics of the J = oo limit 
is not stable against small perturbations. Fig. shows 
the classical spin dynamics of the full model (with T — 
1 and B — 0.1) for a very strong but finite coupling 
J = 100. Here, the motion of the classical spin gets 
very complicated as compared with the highly regular 




FIG. 6: (Color online) Real-time dynamics in the strong-J 
regime: J = 100 and B = 0.1. First panel: Classical spin 
Sit). Second panel: Torque on S{t) due to the exchange 
interaction. Third panel: Torque on 5(t) due to the field 
term. Fourth panel: Angle enclosed by S{t) and {sig)t 


behavior in the weak-J regime (cf. Fig.[^. In particular, 
the z-component of 5(t) is oscillating on nearly the same 
scale as the x and y components. This characteristic 
time scale At « 4 of the oscillation corresponds to a 
frequency oj « 1.5 which differs by more than an order 
of magnitude from both, the Larmor frequency B = 0.1 
and from the exchange-coupling strength J = 100. Note 
that the oscillation of 5, (t) is actually the reason for the 
ambiguity in the determination of the precise reversal 
time and gives rise to the error bars in Fig. [^for strong 
J. 

We attribute the complexity of the dynamics to the 
fact that the torque due to the field term and the torque 
due to the exchange coupling are of comparable magni¬ 
tudes, see second and third panel in Fig.[^ It is interest¬ 
ing that even for strong J, where one would expect 5(t) 
and isig)t to form a tightly bound local spin-zero state, 
there is actually a small deviation from perfect antiparal¬ 
lel alignment, i.e., 7 (t) tt, as can be seen in the fourth 
panel of Fig. This results in a finite z-component of the 
torque on 5(<) which leads to a very fast reversal with 
5z(t) « +0.5 at time t « 2.1. Contrary to the weak- 
coupling limit, however, the z-component of the torque 
changes sign at this point and drives the spin back to the 
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—z-direction. At i « 3.2, however, the z-component of 
S{t) once more reverses its direction. Here, the torque 
due to the exchange coupling vanishes completely as S {t) 
and {si^)t are perfectly antiparallel (see the first zero of 
7 (t) in the fourth panel). The motion continues due to 
the non-zero field-induced torque. This pattern repeats 
several times. S(t) mainly oscillates within a plane in¬ 
cluding and slowly rotating around the z-axis. 

Eventually, there is a perfect relaxation of the classical 
spin for large t but in a very different way as compared to 
the weak-coupling limit. While the deviation from 7 = tt 
is small at any instant of time as for weak J, the most 
apparent difference is perhaps that the new ground state 
is approached with an oscillating behavior of "/(t) around 
7 = TT, i.e., (sio)t may run behind or ahead of S{t) as 
well. 

Let us summarize at this point the main differences 
between the quantum-classical hybrid and the effective 
LLG dynamics of the classical spin: For weak J and B, 
the qualitative behavior, processional motion and relax¬ 
ation, is the same in both approaches. Quantitatively, 
however, the LLG equation is inconsistent with the ob¬ 
served J dependence of the reversal time when assuming 
a oc J^. The B dependence of l/r is linear as expected 
from the LLG approach. For strong J, the spin dynamics 
qualitatively differs from LLG dynamics and gets more 
complicated with a new time scale emerging. Absence of 
complete relaxation, as observed in the strong-H limit, is 
also not accessible to an effective spin-only theory. 


B. Energy dissipation 


To complete the picture of the relaxation dynamics 
of the classical spin, the discussion should also comprise 
the dynamics of the electronic degrees of freedom. The 
spin relaxation must be accompanied by a dissipation of 
energy and spin into the bulk of the electronic system 
since the total energy and the total spin are conserved 
quantities, see Eqs. (12) and (13), while conservation of 


the total particle number, Eq. (11), is trivially ensured by 


the particle-hole symmetric setup considered here where 
the average conduction-electron number at every site is 
time-independent: = 1 - 


The total energy is given by Ftot = {H), see Eq. (13), 
and is a sum over different contributions, Etot = E^{t) 4- 
A'hop(l) + Aint (t), namely the interaction energy with the 
field 


E^{t) = -BS{t ), (16) 

the kinetic (hopping) energy of the conduction-electron 
system 

Ahop(t) = -T^5](cLc,,,)*, (17) 

(b'> ^ 


and the exchange-interaction energy 


Eiretit) = J{sig)tS{t) . 


(18) 




FIG. 7: (Color online) Different contributions to the total 
energy as functions of time for J = 5 and i? = 1. Top panel: 
Eb, energy of the classical spin in the external field. Second 
panel: Ehop/L, kinetic (hopping) energy of the conduction- 
electron system per site, and EtotIL, total energy per site. 
Third panel: Ei^t, exchange-interaction energy. The fourth 
panel shows the time-dependence of the total-energy density 
at different distances from the site io- 


The time dependence of those contributions is shown in 
Fig. for J = 5 and B = 1. 

The top panel of Fig. [7] shows that the system re¬ 
leases the interaction energy | 2 BS'| of the classical spin 
in the external field by aligning the spin to the field di¬ 
rection. In the long-time limit, this energy is stored in 
the conduction-electron system: The average kinetic en¬ 
ergy per site {L = 1001) increases by the same amount 
as shown in the second panel (note the different scales). 
The exchange-interaction energy changes with time but 
is the same for t = 0 and t f 00 (see third panel). The 
total energy is constant (second panel). 

The relaxation of the classical spin in the external field 
B implies an energy flow away from the site ig into the 
bulk of the conduction-electron system such that locally, 
in the vicinity of zq the system is in its new ground state. 
To discuss this energy flow, it is convenient to consider 
the total energy as a lattice sum Etot = ^he 
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FIG. 8: (Color online) Spatiotemporal evolution of the total- 
energy density ei{t), as defined in Eq. (|19[). Calculation for 


J = 5, B = 1. 


total energy “density” defined as 

n.n.(2) 

=—T + Siig{J{sig)t — B)S{t) , 

3 

(19) 

where the sum over j runs over the nearest neighbors of 
site i. The time dependence of the energy density in the 
vicinity of ig and at distances 50 and 100 is shown in the 
fourth panel of Fig. 

At any site in the conduction-electron system, the 
energy density increases from its ground-state value, 
reaches a maximum and eventually relaxes to the en¬ 
ergy density of the new ground state. Since the latter 
is just the ground state with the reversed classical spin, 
the new ground-state energy density is the same as in 
the initial state at t = 0. As can be seen in the fourth 
panel of Fig. there is also a slight spatial oscillation 
of the ground-state energy density which just reflects the 
Friedel oscillations around the impurity at fg. 

Complete relaxation means that the excitation energy 
is completely removed from the vicinity of ig and trans¬ 
ported into the bulk of the system. That this is in fact 
the case can be seen by comparing the energy density at 
different distances from fg. It is also demonstrated by 
Fig. [8] which visualizes the energy-current density which 
symmetrically points away from ig: The total energy of 
the excitation flowing through each pair of sites fg ± Ai is 
constant, i.e., the time-integrated energy flux, f dtei(t) 
is the same for all i. 

As is seen in Fig. (fourth panel) and Fig. there is 
a considerable dispersion of the excitation wave packet 
carrying the energy. For example, at ig -|- 100 it takes 
more than four times longer, as compared to ig -I- 1, un¬ 
til most of the excitation has passed through (note the 
logarithmic time scale). 

The broadening of the wave packet, due to dispersion, 
is asymmetric and bound by an upper limit for the speed 
of the excitation which is roughly set by the F ermi veloc¬ 
ity vf = 2T. This Lieb-Robinson bouncP^*^ determines 
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FIG. 9: (Color online) Spatiotemporal evolntion of the total- 
spin density (si)t (upper panel: ^-component, lower panel 2 - 
component) for J = 5 and B = 0.1. See Fig. |10| for snapshots 
at times indicated by the arrows. 


the “light cone” seen in Fig. 


C. Spin dissipation 

The same upper speed limit, given by the Fermi ve¬ 
locity of the conduction-electron system, is also seen in 
the spatiotemporal evolution of the conduction-electron 
spin density {si)t. This is shown in Fig. for a different 
magnetic field strength B = 0.1 where the classical spin 
dynamics is slower. Apparently, the wave packet of exci¬ 
tations emitted from the impurity not only carries energy 
but also spin. It symmetrically propagates away from Ig 
and, at < « 300, reaches the system boundary where it 
is reflected perfectly. Up to t = 500 there is hardly any 
effect visible in the local observables close to fg that is 
affected by the finite system size. 

Snapshots of the conduction-electron spin dynamics 
are shown in Fig. for the initial state at f = 0 and for 
states at four later times t > 0 which are also indicated 
by the arrows in Fig. At f = 0 the conduction-electron 
system is in its ground state for the given initial direction 
of the classical spin. The latter basically points into the 
—2 direction, apart from a small positive a;-component 
(■d = 7r/50) which is necessary to break the symmetry 
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sites 

FIG. 10: (Color online) Snapshots the of conduction-electron magnetic moments {si}t at different times t as indicated on the 
right and by the corresponding arrows in Fig. Red lines: «-components of {si}t- Blue lines: x components. The profiles are 
perfectly symmetric to the impurity site i = io but displayed up to distances \i — io\ < 100 on the left-hand side and up to the 
system boundary, \i — io| < 500, on the right-hand side. Parameters J = 5, B = 0.1. 


of the problem and to initiate the dynamics. This tiny 
effect will be disregarded in the following. 

From the perspective of the conduction-electron sys¬ 
tem, the interaction term JSsi^ acts as a local external 
magnetic field JS which locally polarizes the conduction 
electrons at iq. Since J is antiferromagnetic, the local 
moment {si^) points into the +z direction. At half-filling, 
the conduct ion-electron system exhibits pronounced an¬ 
tiferromagnetic spin-spin correlations which give rise to 
an antiferromagnetic spin-density wave structure aligned 
to the z axis at t = 0, see first panel of Fig. [Tol 

The total spin Stot = 0 at t = 0, i.e., the classical 
spin S is exactly compensated by the total conduction- 
electron spin (stot) = ground state. 

This can be traced back to the fact that for a I? = 1- 
dimensional tight-binding system with an odd number of 
sites L, with N = L and with a single static magnetic 
impurity, there is exactly one localized state per spin pro¬ 
jection a, irrespective of the strength of the impurity po¬ 
tential (here given by JS = O.SJe^). The number of f 
one-particle eigenstates therefore exceeds the number of 
I states by exactly one. 

Since the energy of the excitation induced by the ex¬ 
ternal field B is completely dissipated into the bulk, the 
state of the conduction-electron system at large t (but 
shorter than t « 500 where finite-size effects appear) 
must locally, close to io, resemble the conduction-electron 
ground state for the reversed spin S — -|-0.5e^. This 
implies that locally all magnetic moments {si)t must re¬ 
verse their direction. In fact, the last panel in Fig. 


(left) shows that the new spin configuration is reached 
for t = 250 at sites with distance \i — io\ ^ 100, see 
dashed line, for example. For later times the spin con¬ 
figuration stays constant (until the wave packet reflected 
from the system boundaries reaches the vicinity of io)- 
The reversal is almost perfect, e.g., {sig)t=o = 0.2649 —>■ 


(sio)t >250 = —0.2645. Deviations of the same order 
of magnitude are also found at larger distances, e.g., 
i = ip —100. We attribute those tiny effects to a weak de¬ 
pendence of the local ground state on the non-equilibrium 
state far from the impurity at t = 250, see right part of 


the last panel in Fig. 10 


The other panels in Fig. demonstrate the mecha¬ 
nism of the spin reversal. At short times (see t = 60, 
second panel) the perturbation of the initial equilibrium 
configuration of the conduction-electron moments is still 
weak. For < = 80 and t = 100 one clearly notices the 
emission of the wave packet starting. Locally, the an¬ 
tiferromagnetic structure is preserved (see left part) but 
superimposed on this, there is an additional spatial struc¬ 
ture of much longer size developing. This finally forms 
the wave packet which is emitted from the central re¬ 
gion. Its spatial extension is about A « 300 as can be 
estimated for t = 250 (last panel on the right) where 
it covers the region 200 < i < 500. The same can be 
read off from the upper part of Fig. Assuming that 
the reversal of each of the conduction-electron moments 
takes about the same time as the reversal of the classi¬ 
cal spin, A is roughly given by the reversal time times 
the Fermi velocity and therefore strongly depends on J 
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and B. For the present case, we have ri « 150/T which 
implies A « 150 x 2 = 300 in rough agreement with the 
data. 

In the course of time, the long-wave length structure 
superimposed on the short-range antiferromagnetic tex¬ 
ture develops a node. This can be seen for t = 100 and 
i « 40 (fourth panel, see dashed line). The node marks 
the spatial border between the new (right of the node, 
closer to io) and the original antiferromagnetic structure 
of the moments and moves away from ig with increasing 
time. 

At a fixed position i, the reversal of the conduction- 
electron moment (sj)t takes place in a similar way as 
the reversal of the classical spin (see both panels in Fig. 
1^ for a fixed i). During the reversal time, its x and y 
components undergo a processional motion while the z 
component changes sign. Note, however, that during the 
reversal | (s^) | gets much larger than its value in the initial 
and in the final equilibrium state. 


V. EFFECTIVE CLASSICAL SPIN DYNAMICS 
A. Perturbation theory 

Eqs. (jrj) and do not form a closed set of equations of 
motion but must be supplemented by the full equation of 
motion ( |10[ ) for the one-particle conduction-electron den¬ 
sity matrix. This implies that the fast electron dynamics 
must be taken into account explicitly even if the spin 
dynamics is much slower. Hence, there is a strong mo¬ 
tivation to integrate out the conduction-electron degrees 
of freedom altogether and to take advantage from a much 
larger time step within a corresponding spin-only time- 
propagation method. Unfortunately, a simple effective 
spin-only action ca n be obtained in the weak-coupling 
(small-J) limit only.E^ES This weak-coupling approxima¬ 
tion is also implicit to all effective spin-only approaches 
that consider the effect of conduction electrons on the 
spin dynamics.^ 

In the weak-J limit the electron degrees of freedom can 
be eliminated in a straightforward way by using standard 
linear-response theoryWe assume that the initial state 
at t = 0 is given by the conduction-electron system in its 
ground state or in thermal equilibrium and an arbitrary 
state of the classical spin. This may be realized formally 
by suddenly switching on the interaction J[t) at time 
t = 0, i.e., J{t) = J0(t) and by switching the local field 
from some initial value Bini at t = 0 to a final value B 
for t > 0. The response of the conduction-electron spin 
at io and time f > 0 {{sig)t = 0 for f = 0) due to the 
time-dependent perturbation J{t)S(t) is 

{s,,)t = J f Sit') (20) 

Jo 

up to linear order in J. Here, the free (J = 0) local 
retarded spin susceptibility of the conduction electrons 


t') is a tensor with elements 


it, t') = -^eit - t')i[si it), 4 (t')]), (21) 

where a,(3 = x,y,z. Using this in Eq. ([^, we get an 
equation of motion for the classical spins only. 


dt 


Sit) = Sit) X B 


-J^S'(t) X f dt'll^'''^^\t-t') ■ Sit') 
Jo 


( 22 ) 


which is correct up to order J^. 

This represents an equation of motion for the classical 
spin only. It has a temporally non-local structure and 
includes an effective interaction of the classical spin at 
time Sit) with the same classical spin at earlier times 
t' < t. In the full quantum-classical theory where the 
electronic degrees of freedom are taken into account ex¬ 
actly, this retarded interaction is mediated by a non¬ 
equilibrium electron dynamics starting at site io and time 
t' and returning back to the same site io at time t > t'. 
Here, for weak J, this is replaced by the equilibrium and 
homogeneous-in-time conduction-electron spin suscepti¬ 
bility — t'). Compared with the results of the 

full quantum-classical theory, we expect that the pertur¬ 
bative spin-only theory breaks down after a propagation 
time t ~ 1/J at the latest. 

Using Wick’s theorem,!^ the spin susceptibility is eas¬ 
ily expressed in terms of the greater and the lesser equi¬ 
librium one-particle Green’s functions, ^^,it,t') = 

-i{ciait)cl^,it')) and G)= (t, f) = i{cl^,it')c^ait)), re¬ 

spectively: 




xlm tr2x2 


a'^G>,(t,t')a'^'GGit',t) 


(23) 


Assuming that the conduction-electron system is charac¬ 
terized by a real, symmetric and spin-independent hop¬ 
ping matrix T^- (as given by the first term of Eq. 

and G'' are unit matrices with respect to the spin 
indices. They are easily expressed as explicit functions 
of T (see Ref.[5Tl for example). With tr (cr“tT“') = 26aa', 
we find 


Tl‘'Z"Ht-t') = Qit-t')S^o,'lm 


1 -|- 






(24) 


* 0*0 


^O'lQ 


for a conduction-electron system at inverse temperature 
P and chemical potential /i. 

Using Eq. (24) we have computed the spin suscepti¬ 


bility for the ground state (/3 = oo) of the conduction- 
electron system. This fixes the kernel in the integro- 
differential equation Eq. \22\ which is solved numeri¬ 
cally by standard techniques.^ We again consider the 
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FIG. 11: (Color online) Components Sx{t) and Sz{t) of the 
classical spin after a sudden switch of the field from x to 
2 ; direction at t = 0. Calculations for B = 1 and J = 1. 
Solid lines: results of the linear-response dynamics, Eq. (22l. 
Dashed lines: results of the exact quantum-classical dynamics 
for L = 1001. 


system displayed in Fig. with a single classical spin 
coupled via J to the central site jq of a chain consisting 
of L = 1001 sites. Finite-size artifacts do not show up 
before tmax = 500. 

Figs, pd] and 12 show the resulting linear-response dy¬ 
namics of the classical spin after preparing the initial 
state of the system with the classical spin pointing into 
+x direction while B = (0,0, B). The external magnetic 
field induces a precessional motion of the classical spin: 
there is a rapid oscillation of its x component (and of its 
y component, not shown) with frequency to ~ B (blue 
lines). Damping is induced by dissipation of energy and 
spin: for large times, the 2 component aligns to the ex¬ 
ternal field (red lines). 

For weak coupling, up to J = 1 (Fig.[^, there is an al¬ 
most perfect agreement between the results of the exact 
quantum-classical dynamics (full lines) and the linear- 
response theory (dashed lines) up to the maximum prop¬ 
agation time tmax = 500. We note that, compared to 
the full theory, there is a tiny deviation of the linear- 
response result for the z component of S(t) visible in 
Fig. 11 for times t > 10. Hence, on this level of accuracy, 
ti « 10 sets the time scale up to which the linear-response 
theory is valid. This may appear surprising as this im¬ 
plies J = 10 for the “small” dimensionless parameter 
of the perturbation theory. One has to keep in mind, 
however, that even if the perturbation is “strong”, its 
effects can be rather moderate since only non-adiabatic 
terms ~ S{t) x S{t') contribute in Eq. (221. 

For J = 3, see Fig. damping of the classical spin 
sets in much earlier. Visible deviations of the linear- 
response theory from the full dynamics already appear 
on a time scale that is almost two orders of magnitude 
smaller as compared to the case J = 1. A simple rea¬ 
soning based on the argument that the dimensionless 
expansion parameter is ti J fails as this disregards the 
strong enhancement of retardation effects with increas¬ 


FIG. 12: (Color online) The same as Fig. 11 but for a different 
exchange coupling constant J = 3. 


ing J, which have been discussed in Sec. |1V A[ These 
effects make the perturbation much more effective, i.e., 
lead to a torque, which is exerted by the conduction elec¬ 
trons on the classical spin, growing stronger than linear 
in J. 

We conclude that linear-response theory is highly at¬ 
tractive formally as it provides a tractable spin-only effec¬ 
tive theory. On the other hand, substantial discrepancies 
compared to the full (non-perturbative) theory show up 
as soon as damping effects become stronger. Note that, 
with increasing time, these deviations must diminish and 
disappear eventually since both, the full and the effec¬ 
tive theory, predict a fully relaxed spin state for t —00 
- see lower panel of Fig. for example. At least for 
simple systems with a single classical spin, as considered 
here, this implies that the effective theory provides qual¬ 
itatively reasonable results. 


B. Landau-Lifschitz-Gilbert equation 


To derive the LLG equation , the linear-response the¬ 
ory mus t be further simplifiedi ^^l^^l As is obvious from 
Eq. (22), the spin-susceptibility —t') can be inter¬ 


preted as an effective retarded self-interaction of the spin. 
We assume that the electron dynamics is much faster 
than the spin dynamics. On the time scale of the spin 
dynamics, the self-interaction then takes place almost in¬ 
stantaneously, i.e., the memory kernel — t') = 

— t') in Eq. (22) is peaked at t' « t. We can 
therefore approximate S{t') « S{t) + {t' — t)S{t) under 


the integral in Eq. (22). This immediately yields: 


dSjt) 

dt 


= S{t) xB + a{t)S{t) X 


dSjt) 

dt 


(25) 


where 


a{t) = 



(26) 
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after substituting t' ^ t = t — t'. 

Eq. (25) takes the form of the standard LEG equation 
for a single classical spin [cf. Eq. §] if a{t) is replaced 
by 


noo 

lim a{t) = J^ . 


(27) 


Note that this is a necessary step to arrive at a constant 
damping parameter which can again be justified by not¬ 
ing that nb®*)(t) is peaked at t = 0. 

Before proceeding, let us stress that Eq. ( p^ is, or 
is equivalent to, the standard expression used for com¬ 
puting the Gilbert damping constant in various studies: 
After Fourier transformation, 


nb®‘)(a;) =y' , 


(28) 


one ends up with 

a = -ij2^nb®‘)(w) 
aw 


= J 2 ^ 1 mnb®‘)(w) 

cj^O OUJ 


(29) 

which has also been derived, e.g., in Refs. I14I54I in dif¬ 
ferent contexts. The freq uency-dependent spin correla¬ 
tion nb®*)(cLi) in Eq. (29) can be obtained explicitly as 
the Fourier transform of nb®*)(t) given by Eq. (24). A 
straightforward calculation yields: 

d^ df(w) ^ (- 30 ) 


^ t 2 

a=-J 


dw 


where f{w) = l/(exp(/3tLi) -|-1) is the Fermi function, and 


M I I I I I I I I M l| I I I I I I I 11 



time t 


FIG. 13: (Color online) Local retarded spin susceptibility of 
the conduction electrons 11*''®*^ (t) (red line) and 
(blue) as obtained from Eq. ( |33[ ) for a one-dimensional 
conduction-electron system with L = 10000 sites and peri¬ 
odic boundary conditions. 


C. Ill-defined Gilbert damping 


The above discussion shows that Eq. (27) represents 


the fundamental definition of the Gilbert damping con¬ 
stant a and that the limit f —>■ 00 is crucial to recover the 
LEG equation in its standard form. The existence of the 
long-time limit, however, decisively depends on the long¬ 
time behavior of the retard ed s pin-correlation function 
n(ret)(^)^ Starting from Eq. (24), this is easily computed 


as 


Aioc(w)=L ^ y^(5(a; -b /r - e{k)) 


(31) 


for L —>■ 00 is the local one-particle spectral function. 
Here we have assumed periodic boundary conditions, i.e., 
spatial homogeneity, such that the hopping matrix T is 
diagonalized by Fourier transformation: 




kj 


(32) 


with Uik = exp{ikRi). The eigenvalues of T 

are given by the tight-binding dispersion e{k) of the 
conduction-electron Bloch band. 

Within our simple tight-binding model for the 
conduction-electron system, Eqs. (29) and (30) are equiv¬ 


alent with Kambersky’s breathing Fermi-surface the¬ 
ory, related torque-correlation models and scattering 
theory and have frequently been used for ab initio as 
well as model computations of the Gilbert damping 

constant.™l2ni2lIlMI] 


Eet us remark that Eq. (30) demonstrates that a < 0. 


This results from the convention for the coupling of the 
magnetic held to the spin, namely H = H{B = 0) — BS 
[see Eq. i)]. which has been adopted here. As a conse¬ 
quence, the precession of S{t) around B is described by 
a le/thand helix, S = S{t) x B, and thus a must be 
negative to describe damping. 


n(-‘)(f) = 0(t)lm_^ 




Je{p)t 


L2 I g-^(e(fc)-M) e/3(e(p)-Ai) -|- 1 ' 
k,p 

(33) 

Fig. gives an example for the time-dependence of 
n(ret) ^ rpj^g calculations have been done at half-hlling, 
/? = 00 for L = 10^ sites. We note that the susceptibility 
is in fact peaked at t = 0. For long times, it oscillates 
with frequency wn = 4 and decays as 1/t. This is an im¬ 
portant observation as it implies that the limit in Eq. (27) 


does not exist and that, therefore, the damping constant 
a is ill-dehned. 

To analyze the physical origin of the divergent integral, 
we rewrite the spin susceptibility in Eq. ([M|) as 


nM)(i) = 


^loc 




(34) 


Its long-time behavior is governed by the long-time be¬ 
havior of the Fourier transform 


^(occ,unocc)^^^^ _ / ^ ^ (occ,unocc) 


^loc 


\t)= / 


loc 


iw) (35) 


of the occupied, a|°“^(w) = /(a;)Aioc(a;), and of the un¬ 
occupied, (w) = f{w)A\oc{w), part of the spectral 

density, see Eq. (31). 


For functions with smooth w dependence, the Fourier 
transform generically drops to zero exponentially fast if 



































































13 


i > oo. A power-law decay, however, is obtained if there 
are singularities of We can distinguish 

between van Hove singularities, which are, e.g., of the 
form oc Q{uj — uJo){uj — ujq)'^ (with k > — 1 ), and the step¬ 
like singularity oc 0(a; — wq) (i-e. k = 0 ), arising in the 
zero-temperature limit at wq = 0 due to the Fermi func¬ 
tion. Generally, a singularity of order k gives rise to the 
asymptotic behavior apart from 

a purely oscillatory factor For the present case, 

the van Hove singularities of at iwo = 2 

explain, via Eq. (34) the oscillation of with fre¬ 

quency wn = 2a;o = 4. 

Generally, the location of the van Hove singularity 
on the frequency axis, i.e. wq, determines the oscilla¬ 
tion period while the decay of is governed by 

the strength of the singularity. Gonsider, as an exam¬ 
ple, the zero-temperature case and assume that there are 
no van Hove singularities. The sharp Fermi edge im¬ 
plies oc The 

Gilbert-damping constant is well defined in this case. 

The strength of van Hove singularities depends on 
the lattice dimension For a one-dimensional lat¬ 

tice, we have van Hove singularities with k = — 1/2, and 
thus cx consistent with Fig. (13). Here, 

the strong van Hove singularity dominates the long-time 
asymptotic behavior as compared to the weaker Fermi- 
edge singularity. For D — 3, we have k = 1/2 and 
n(ret)(^) Qj. ^-3 if ^ < (X) while for /3 = C !0 the Fermi- 
edge dominates and oc The D = 2 case is 

more complicated: The logarithmic van Hove singularity 
cx ln|w| leads to n(’'®*)(t) oc This, however, applies 
to cases off half-filling only. At half-filling the van Hove 
and the Fermi-edge singularity combine to a singularity 
cx ©(w) In |a;| which gives n("'®*)(t) cx ln^(t)/t^. For finite 
temperatures, we again have cx 

The existence of the integral Eq. (271 depends on the 
t —)■ oo behavior and either requires a decay as H^’'®*^ (t) cx 
or faster, or an asymptotic form n(‘'®*^(t) oc 
with an oscillating factor resulting from a non-zero po¬ 
sition c^o 7 ^ 0 of the van Hove singularity. For the one¬ 
dimensional case, we conclude that the LEG equation 
(with a time-independent damping constant) is based on 
an ill-defined concept. Also, the derivation of Eqs. (29) 
and (30) is invalid in this case as the ui derivative and the 
t integral do not commute. This conclusion might change 
for the case of interacting conduction electrons. Here one 
would expect a regularization of van Hove singularities 
due to a finite imaginary part of the conduction-electron 
self-energy. 


VI. CONCLUSIONS 

Hybrid systems consisting of classical spins coupled 
to a bath of non-interacting conduction electrons rep¬ 
resent a class of model systems with a non-trivial real¬ 
time dynamics which is numerically accessible on long 


time scales. Here we have considered the simplest vari¬ 
ant of this class, the Kondo-impurity model with a clas¬ 
sical spin, and studied the relaxation dynamics of the 
spin in an external magnetic field. As a fundamental 
model this is interesting of its own but also makes con¬ 
tact with different fields, e.g., atomistic spin dynamics in 
magnetic samples, spin relaxation in spintronics devices, 
femto-second dynamics of highly excited electron systems 
where local magnetic moments are formed due to electron 
correlations, and artificial Kondo systems simulated with 
ultracold atoms in optical lattices. 

We have compared the coupled spin and electron dy¬ 
namics with the predictions of the widely used Landau- 
Lifshitz-Gilbert equation which is supposed to cover the 
regime of weak local exchange J and slow spin dynamics. 
For the studied setup, the LEG equation predicts a rather 
regular time evolution characterized by spin precession, 
spin relaxation and eventually reversal of the spin on a 
time scale r depending on J (and the field strength B). 
We have demonstrated that this type of dynamics can be 
recovered and understood on a microscopic level in the 
more fundamental quantum-classical Kondo model. It is 
traced back to a non-adiabatic dynamics of the electron 
degrees of freedom and the feedback of the electronic sub¬ 
system on the spin. It turns out that the spin dynamics 
is essentially a consequence of the retarded effect of the 
local exchange. Namely, the classical spin can be seen as 
a perturbation exciting the conduction-electron system 
locally. This electronic excitation propagates and feeds 
back to the classical spin, but at a later time, and thereby 
induces a spin torque. 

We found that this mechanism drives the relaxation of 
the system to its local ground state irrespective of the 
strength of the local exchange J. As the microscopic dy¬ 
namics is fully conserving, the energy and spin of the 
initial excitation which is locally stored in the vicinity of 
the classical spin, must be dissipated into the bulk of the 
system in the course of time. This dissipation could be 
uncovered by studying the relaxation process from the 
perspective of the electron degrees of freedom. Dissipa¬ 
tion of energy and spin takes place through the emission 
of a dispersive spin-polarized wave packet propagating 
through the lattice with the Fermi velocity. In this pro¬ 
cess the local conduction-electron magnetic moment at 
any given distance to the impurity undergoes a reversal, 
characterized by precession and relaxation, similar to the 
motion of the classical spin. 

The dynamics of the classical spin can be qualitatively 
very different from the predictions of the EEG equation 
for strong J. In this regime we found a complex mo¬ 
tion characterized by oscillations of the angle between 
the classical spin S{t) and the local conduction-electron 
magnetic moment at the impurity site (sig) around the 
adiabatic value 7 = tt which takes place on an emergent 
new time scale. 

In the weak-J limit, the classical spin dynamics is qual¬ 
itatively predicted correctly by the EEG equation. At 
least partially, however, this must be attributed to the 
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fact that the LLG approach, by construction, recovers 
the correct final state where the spin is parallel to the 
field. In fact, quantitative deviations are found during 
the relaxation process. The LLG approach is based on 
first-order perturbation theory in J and on the additional 
assumption that the classical spin is slow. To pinpoint 
the source of the deviations, we have numerically solved 
the integro-differential equation that is obtained in first- 
order-in-J perturbation theory and compared with the 
full hybrid dynamics. The deviations of the perturbative 
approach from the exact dynamics are found to gradually 
increase with the propagation time (until the proximity 
to the final state enforces the correct long-time asymp¬ 
totics). This is the expected result as the dimensionless 
small parameter is Jt. However, with increasing J the 
time scale on which perturbation theory is reliable de¬ 
creases much stronger than 1/J due to a strong enhance¬ 
ment of retardation effects which make the perturbation 
more effective and produce a stronger torque. 

Generally, the perturbation can be rather ineffective in 
the sense that it produces a torque oc S{t) x S{t') which 
is very weak if the process is nearly adiabatic. This ex¬ 
plains that first-order perturbation theory and the LLG 
equation is applicable at all for couplings of the order of 
hopping J ~ T. For the present study this can also be 
seen as a fortunate circumstance since the regime of very 
weak couplings J ^ T is not accessible numerically. In 
this case the spin-reversal time scale gets so large that 
the propagation of excitations in the conduction-electron 
subsystem would by affected by backscattering from the 
edges of the system which necessarily must be assumed 
as finite for the numerical treatment. 

For the one-dimensional lattice studied here, a di¬ 
rect comparison between LLG equation and the exact 
quantum-classical theory is not meaningful as the damp¬ 
ing constant a is ill-defined in this case. We could argue 
that the problem results from the strength of the van 


Hove singularities in the conduction-electron density of 
states which dictates the long-time behavior of the mem¬ 
ory kernel of the integro-differential equation which is 
given by the equilibrium spin susceptibility. As the type 
of the van Hove singularity is characteristic for all sys¬ 
tems of a given dimension, we can generally conclude that 
the LLG approach reduces to a purely phenomenological 
scheme in the one-dimensional case. However, it is an 
open question, which will be interesting to tackle in the 
future, if this conclusion is still valid for systems where 
the Coulomb interaction among the conduction electrons 
is taken into account additionally. 

There are more interesting lines of research which are 
based on the present work and could be pursued in the 
future. Those include systems with more than a single 
spin where, e.g., the effects of a time-dependent and 
retarded RKKY interaction can be studied additionally. 
We are also working on a tractable extension of the the¬ 
ory to account for longitudinal fluctuations of the spins 
to include time-dependent Hondo screening, and the 
competition with RKKY coupling, on a time-dependent 
mean-field level. Finally, lattice rather than impurity 
variants of the quantum-classical hybrid model are 
highly interesting to address the time-dependent phase 
transitions. 
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